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Q\ ', Introduction 

A quantum many particle system on the line with an interaction potential U(x) 
is defined by the Hamiltonian 

> ■ 

■ N d 2 

(1) H = J2d^ + K 2 U(xi-Xj), KeR. 

(N ' j=l 3 Ki<j<N 

(N ; 

Such systems were first considered by Calogero for U (x) = x~ 2 [Ca] and Sutherland 
j> . for U(x) = sinh(a;) -2 [Su]. 

Any differential operator commuting with H is called a quantum integral of the 
\ system. One says that two quantum integrals are in involution if they commute as 
^-i ■ differential operators. 

^ ■ The system defined by H is called completely integrable if it has A quantum 

integrals in involution, L\, L^, which are algebraically independent: every poly- 
nomial P G M[zi, Zn] such that P(Li, L^) = is identically zero. 

For A > 2, the system defined by H is not always completely integrable. How- 
ever, for special choices of U(x) it is known that it is completely integrable for any 

^ ■ K. For example, this is proved for U(x) = x~ 2 (rational case), U(x) = sinh(x) -2 

> . (trigonometric case), and U(x) = p(x\r), where p is the Weierstrass elliptic func- 
k> ! tion: 



(2) p(x\t) = ^+ E 

(m,n)€Z 2 \{(0,0)} 

(elliptic case) (see [OP] and references therein) . These results come from an explicit 
construction of quantum integrals via the Lax matrix. 

In this paper we propose a new construction of the quantum integrals for H in 
the trigonometric and elliptic cases using representation theory of Lie algebras, and 
give a new proof of the complete integrability theorem based on this construction. 
In the trigonometric case, we construct a homomorphism from the center of the 
universal enveloping algebra of qI n to the algebra of differential operators in N 
variables. We show that the image of the second-order Casimir under this homo- 
morphism is a multiple of the Sutherland operator, which implies that the images of 
higher Casimirs are its quantum integrals and thus proves the integrability. In the 
elliptic case, we apply a similar method to the center of the (completed) universal 
enveloping algebra of the affine gl N at the critical level k = —A, and thus produce 
quantum integrals of the Hamiltonian (1) with the elliptic potential U = Kp. 

If U(x) is one of the above three choices, we call H a Calogero-Sutherland (CS) 
operator. If H is a CS operator then for any set of complex numbers Aj the system 
of differential equations 
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is consistent and has AH linearly independent solutions. This system is called the 
eigenvalue problem for H. The solutions of (3) for the rational and trigonometric 
cases (unlike the elliptic case) have been studied in many papers [*] and are rather 
well understood. We will give a new construction of these eigenfunctions for the 
trigonometric case as normalized traces of intertwiners between representations of 
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1. Diagonalization of the Sutherland operator 



The Sutherland operator is the following differential operator: 

N f)2 K 

(L1) E sinV{xi . Xj y K ^ 

j = l 3 l<i<j<N v 1 31 

The quantum integrals for this operator and their eigenfunctions (zonal spherical 
functions) are known [OP,HO]. In this section we will give a new description of the 
quantum integrals and eigenfunctions using representatin theory of gl N . 

Let qI n denote the Lie algebra of complex N x N matrices. If A e gl N , we will 
denote by Aij the entry at the intersection of the i-th row and j-th column of A. 
Also, Eij will denote the elementary matrices: (Eij) mn = 5i m 5j n . For brevity we 
denote En by hi. 

Let A = (Ai,...,Ajv) € C N , and let M\ be the Verma module over gl N with 
highest weight A, i.e. the module generated by a highest weight vector v\ satisfying 
the defining relation 

N 

(1.2) Av x = J2 X 3 A n v \i 

i=i 

whenever A is an upper triangular matrix. Let denote the restricted dual 
module to M\ - the direct sum of dual spaces to the weight subspaces in M\ with 
the action of gl N defined by duality. 

Let fx e C. Define a module V M over gl N as follows. As a vector space, V M 
is the space of functions of the form /(x) = (ILjli •••> g" 1 ); where p £ 

Cfyf 1 , . . . , y^-i] is a Laurent polynomial, and the action <fi of gi N in is described 
by the formula 4>(Eij) = — H&ij- 

It is obvious that the set of weights of V M is {(Ai, Xn) E Z n = 0}, and 

every weight occurs with multiplicity 1. 

Consider the completed tensor product M\®V^, = Homc(M^,V^). It has a 
natural structure of a gl^-module. 

Proposition 1.1. If M\ is irreducible then there exists a unique up to a factor 
nonzero intertwining operator $a : M\ — > M\®V^. 

Proof. We need to prove that the module Homc(M^, V M ) contains a singular vec- 



3 



unique. Since w\ is invariant under the Lie subalgebra n + of strictly upper trian- 
gular matrices, it is uniquely determined by its value on the lowest weight vector 
vl of the module M A . Because M\ is irreducible, = U(n + )v^, and therefore 
w\(v^) can be any zero weight vector in V M . Since the zero weight subspace of V M 
is one dimensional, the vector w\ exists and is unique up to a factor, Q.E.D. ■ 

We fix the normalization of the operator $a under which < u A , <&v\ >= Yl£j e 
Let p = (pi, pat), pj = — j. Consider the function of N variables 

( } A( "" N) = TrU, z^Y' 

This function takes values in and has weight zero, hence 

(i.4) *A=^n^ 

where ifj\ is a complex-valued function. 

Let us introduce new coordinates Xi such that Zi = e 2x \ and regard the function 
V>a as a function of (x\, ...,xn) 



If a, 6 G C^, let < a, 6 >= Yl!j=i a 3^j^ an< l a2 = < a ' a >• Also, for brevity 



we 



use the notation Z = z 1 1 . . . z N N , and write Tr A for Tr 

Theorem 1.2. /f K m formula (1.1) equals —2p(p + 1) t/ien £/ie function tfj\ is 
an eigenfunction of the differential operator H defined by (1.1) with the eigenvalue 
4(A + p) 2 . 

Proof. The Casimir element 

N 

(1.5) C 2 = J] EyEji 

acts in Ma by multiplication by < A, A + 2p >. Therefore, 

(1.6) TV ^ A ^f ) =< A, A + 2p > * A ( Zl , z„). 

On the other hand, the left hand side of (1.6) can be written in the form ^fj=i *a > 
where 



(1.7) * 



, , = Tr A ($ A £ tJ £ Jt Z) 
Tr- P (^) 



Let us express the terms in terms of \I/a- First of all, consider the function 
F(z±, zn) = Tr- p (Z). This function is the character of the module M_ p , there- 
fore 



^-i y gj 



(n,=i^) 2 



riR , z i dF_(N-l v ^ 
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We consider two cases. 

1) i = j. Then, using (1.8), we obtain 



If d x ~ 



2 

z. 



9 I p 



2) i j. Then we have 



¥1 = F^TrxiQxEijEjiZ) = 
F^TixiEijQxEjiZ) + F- 1 E ij Tr x ($ x E ji Z) = 
F^Trxi^xEjiZEij) + F^E^Ttx^xE^Z) = 

^-F^Ttx^xEj^Z) +F- 1 E tJ Trx($xE Jl Z) = 



-F^Trx^xEijEjiZ) + ^F" 1 Tr A ($ A (£; jJ - E U )Z) 

+F- 1 Ei J Trx($xE Jl Z) = 
(1.10) |^ + F~^{z^ - Zi ^j (F* A ) + F-'E^Trx^xE^Z). 

Formula (1.10) is a linear equation on ~^>\ J . Solving this equation, we obtain 

It remains to compute Tr x($xEjiZ). We have 

Tr A ($ A £^Z) = TixiEj&xZ) + EjiTrxi^xZ) = Tix(Q\ZEji) + FE^ X = 

(1.12) ^Trxi^xEjiZ) + FEjiVx, 

Zi 

which implies 

(1.13) Trxi^xEjiZ) = -^—FE^x- 

Zi Zj 

Combining (1.11) and (1.13), we deduce 

It is easy to see that the operator EijEji acts in the zero weight subspace of by 
multiplication by + 1). Therefore, (1.14) can be rewritten as follows: 

zj - Zi\ J dzj dzi J 



(1-15) + (zj^- - Zi ^-) logF - pi* + 1W 



ZiZj 
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Now, summing up equations (1.9) for all % and equations (1.14) for all i ^ j, and 
using (1.8) and the identity 

("«) E ^ - = Ef^r - (" - 

we obtain 



AT / o \ 2 



Tr A ($ A C 2 Z) ^E/ 9 



(1.17) 



Tr_ p (Z) t[\ dz 



2 
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From formula (1.8) it follows that 

(1.18) F- 1 {z 1 ,...,z N )=det{z- p >). 
Using this identity, we find 

(1.19) E(^)F-'=<„,P>F- 

Substituting (1.19) into (1.17), we get 
(1.20) 

Rewriting equation (1.20) in the new variables Xi such that Zi = e 2x \ we obtain 
(1.21) 



N 



Tr x (^C 2 Z) = 1 ^ ^ + 1) T 1 

Tr_ p (Z) 4^ Ox 2 ^ ri P ^ a 2 Z^sinhS^-^.) 

Comparing (1.6) and (1.21), we finally get 



N f) 2 I 1 
"2) Z^-2m(m + 1)E ^4(A + P ) 



(1.22) > _ - 2//(// + 1) ^ —77277 77T ^ = 4(A + p)^A- 

The theorem is proved. ■ 



2. Quantum integrals of the Sutherland operator 

Let us now find the quantum integrals of the Sutherland operator. For this 
purpose we will use the higher Casimir elements of the Lie algebra qI n . 

Let Y be any element of the universal enveloping algebra U(gl N ). Then we can 
consider the following V^-valued function: 

(2.i) * l( y|, a) -^-.(^f'-*) 



Tr| M _„ (*?'... 4~) 
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Proposition 2.1. There exists a differential operator Cy in z\, ...,zn whose coef- 
ficients are U(gl N ) -valued functions, such that 

(2.2) *x(Y\ Zl , z N ) = C Y ^x{z u z N ). 

Proof. We say that Y is of order < m if it is a sum of monomials including < m 
factors from g{ N . The proof is by induction in the order of Y. For order the 
statement is obvious. If Y is of a positive order m then we may assume that Y is a 
monomial of the form Y = XEij. Since X is a monomial of a lower order than Y, 
we may assume that the operator Cx is already defined. 
Again, we consider two cases. If i = j then 

(2.3) *x(Y\z u z N ) = F~ l Zi— — (Fty\(X\zi, z N )), 

OZi 

so we can set Cy = F -1 o Zi -^- o F o £ X - 

If % 7^ j then let U = [X,Eij]. Clearly, the order of U is < m — 1, so we can 
assume that Ljj has been defined. Then we have 

* A (y) = ^x(EijX) + * X (U) = E i:i * x (X) + -*x(Y) + £u*x = 

z j 

(2.4) (E ij C x + C u )^x+-^x(Y), 
which implies 

(2.5) #aOO = -^—(EtjCx + £u)®x. 

Zj Zi 

This proves the step of induction Q.E.D. ■ 

Let Y E U(qI n ) be of weight zero, i.e. such that = for all i. Then 

^\(Y\zi, z n) = ip\(Y\x\, xn) YliLi ^f' where ip is a complex-valued function. 

Proposition 2.2. If Y is of weight zero then there exists a unique scalar-valued 
differential operator Dy in the variables xi,...,xn whose coefficients are rational 
functions of e 2x % 1 < % < N such that 

(2.6) il) X {Y\x u ...,x N ) = Dyipxix!, ...,x N ) 
for a generic X. 

Proof. If Y is of weight then Cy preserves weight, so its coefficients map the zero 
weight subspace in to itself. Since the zero weight subspace is one-dimensional, 
the coefficients of Cy become scalars when restricted to this space. In this way the 
operator Cy turns into a scalar differential operator Dy which satisfies (2.6). 
To establish the uniqueness, it is enough to prove the following: 
Lemma 2.3. Any differential operator D annihilating for a generic A whose 
coefficients are rational functions in Zi = e 2xi , 1 < i < N must be identically equal 
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Proof. Since the functions are homogeneous in z\, zn, any homogeneous 
component of D also annihilates them. Therefore, we may assume that D is homo- 
geneous. We can also assume that D is of degree (otherwise we can multiply it 
by a suitable power of z±) and that its coefficients are polynomials in yi = Zi/zi+i 
(otherwise we can multiply D by the common denominator of its coefficients) . Let 
us write D in the form 

(2-7) D=J2y?---y r N-iDr 1 ,...,r N - 1 , 

rj >0 

where = Pn rw .i(^, q§^, where Pn r w _i are polynomials. Let 

k be the smallest integer for which there exist r±, rjy-i with J^rj = k such that 
D ri ,...,r N -i 7^ 0. Since ipl e n e 2xi ^+f^C[[y u y N -i]], for this set of Tj we have 

(2.8) J Dr 1 ,...,r JV _ 1 e <2x ' A+p> = 0, x = On, Xjv). 

This implies that p J . li ... /J ,_ 1 (2(A + p)) =0 for a generic A, which means that 
Pri,...,rjv-i i s identically zero - a contradiction. This completes the proof of the 
lemma and the proposition. ■ 

It is well known that the center of U(gl N ) is freely generated by the Casimir 
elements 

N 

(2.9) C m = E hh ■ ■ ■ E jm _ ljm E jmjl , 1 < m < N. 

jl,---,jm = l 

Define the differential operators 

(2.10) L J= D Cj , l<j<N 
Proposition 2.4. 

(i) Li = iY:?=i&;> L * = T+<P>P>- 

(ii) For any 1 < j < N the function ifj\ is an eigenfunction of Lj with an 
eigenvalue Pj(X + p), where pj is a symmetric polynomial of degree j . 

(Hi) [Li, Lj] = for all 

(iv) The symbol of L m is 2~ m J2iLi 

(v) The operators Lj and the polynomials pj are algebraically independent. 

(vi) The operators Li are invariant under permutations of x\, ...,xn- 

Proof, (i) The formula for L\ is obvious. The formula for L2 follows from Theorem 
1.2. 

(ii) The element Ci acts in M\ by multiplication by a scalar polynomially de- 
pendent on A. This scalar is a symmetric polynomial in A + p, since the action 
of Ci is the same in M\ and M <t (a+ p )- p for any permutation a (proof: the first 
module contains the second one when A is a dominant integral weight, and any two 
polynomials coinciding on dominant integral weights coincide identically). Denote 
this symmetric polynomial by pi. Then (2.10) implies that L^x = Pi{\ + p)i^\- 

(iii) Statement (ii) implies that [Li,Lj]i/j\ = for a generic A. According to 
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(iv) It is easy to show that 

N 

(2.11) C m = K + XijEij + Y, 

where Xij, Y are elements of U(qI n ) of order < m — 1. This fact combined with 
the proof of Proposition 2.1 shows that the symbol of L m will be the same as the 
symbol of D^n hm (since the rest of the terms in (2.11) will give lower order 

i=i j 

contributions), i.e. it will equal 2 _m X^=i ~§x^- 

(v) Since the symbols of Lj are algebraically independent (they are elementary 
symmetric functions), so are Lj themselves. Therefore, so are py. if there existed a 
nonzero polynomial Q such that Q(pi, ■■■,Pn) is identically zero, then Q(Li, Ln) 
would annihilate tf)\ for a generic A, which, by Lemma 2.3, would mean that this 
operator is zero, so that from the algebraic independence of Lj one gets Q = 0. 

(vi) Let a be any transposition. Then for any Y one has C, aYa -\ = aC Y o~~ l , 
where C Y is defined by the formula: 

(2.12) (C a f) (a^i), ...,x a(N) ) =£f(x a ( 1 ),...,x a ( N )). 

Therefore, D aYa -\ = D Y . On the other hand, o"C m o" _1 = C m , which implies that 
L m = L m . ■ 

We can now consider the eigenvalue problem (3) associated with the collection 
of operators {Lj}. Recall that this problem is to find a basis of the AH-dimensional 
space of solutions of the system of differential equations Liip = A^, 1 < i < N, 
where Aj are some given complex numbers. 

Proposition 2.5. Let X be a solution of the system of algebraic equations 

(2.13) p t (\ + p)=A l , l<i<N. 

Suppose that Xi + Pi — Xj — Pj is not an integer for any i ^ j. Then the functions 
{il) a (\ +p )-p, o G S n}, (where Sn is the symmetric group) form a basis in the space 
of solutions of system (3). 

Proof. First of all, all the functions {ip a ^\ + p)- p } are defined since all the modules 
M (T (A+p)-p are irreducible. Next, the weights a(X + p) — p are different for different 
a, hence the functions {'0 cr (A+p)_ p } are linearly independent (this follows from the 

asymptotics ip\ ~ YlZj J+P \ Zi/zi+i — > 0, 1 < i < N — 1). And finally, all the 
functions {i/v(A+p)-p} are solutions of (3): ip\ is a solution of (3) by the definition 
of A, and the rest are solutions because the operators Li and the polynomials pi 
are symmetric and therefore X a = a(X + p) — p satisfies (2.13) for any a G Sn- B 

Remark. In 1976 I.Frenkel computed the Laplace's operator in the space of 
functions on (SU(N) x SU(N))/ SU{N)di ag with values in a symmetric power of the 
fundamental representation of SU (N) equivariant with respect to the left diagonal 
action of SU (N) , and found that it was equal to the Sutherland operator (with sin 
instead of sinh) and K = — 2m(m+l), where m is an integer (I.Frenkel, unpublished 

wnrkl Prnnnsitinn 9, ¥\ ran Vip rpffarrlprl as a crpnprali'Zfltinn nf this rpsnlt 
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3. Affine Lie algebras, vertex operators, 
correlation functions, and regular expressions 

Recall the definition of affine Lie algebras. 

The affine Lie algebra gl N is defined as follows: as a vector space 
gl N = gl N <S> C[t, © Cc, and the commutator is given by 

(3.1) [A(t) + ac, B{t) + [3c] = [A(t),B(t)} + Res t=0 (* _1 Tr(A'(*)B(*))), 

where A(t), B(t) are matrix- valued trigonometric polynomials. 

The principal gradation on gl N assigns degree j — i + Nm to the element Eijt m , 
and degree zero to c. The degree of an element a is denoted by deg(a). 

Define the Lie subalgebra go in gl N to be the Lie algebra of all elements A(t) + 
ac G gl N such that Tr (A(t)) is a constant (does not depend on t). This algebra is 
isomorphic to the direct sum of the affine Lie algebra sf/v and the one-dimensional 
abelian Lie algebra. It inherits the grading from q{ n . We denote by g^ the subal- 
gebras in go spanned by the elements of positive (respectively negative) degree. 

We will work with two kinds of modules over go - Verma modules and evaluation 
modules. 

Let A G and V be a go-module. We say that a nonzero vector v G V is of 
weight A if hiV = \v. We say that a nonzero element a G U(g ) is of weight A if 
[hi, a] = \a. 

Let A G and k G C. The Verma module over g is defined by a single 

generator and the relations: 

N 

(3.2) (A(t) + ac)v x , k = (^ A jj(°) X j + ^k)v x , k , A(t) + ac G g 

i=i 

if A(t) is regular at t = 0, and A(0) is upper triangular. This module is a highest 
weight module which is irreducible for generic A and k. 

The principal gradation operator d on M\^ k is defined by the relations: 

(3.3) dv\ : k = 0, [<9, a] = na if a G go, deg(a) = n. 

A vector w G is called homogeneous if it is an eigenvector of d. We say that 

a vector w G M\,k is of degree n (deg(w) = n) if dw = nw. 

It is easy to check that the only vectors of degree in are multiples of v\^, 
and the rest of homogeneous vectors have negative degrees. Also, the space of vec- 
tors of a fixed degree n is finite dimensional. This space is called the homogeneous 
subspace of degree n, and denoted by M\^[n]. 

We will need to work with the completion of the Verma module with 
respect to the principal gradation. is defined as the space of infinite sums of 

the form J2 n<0 w n , where w n is a vector of degree n. Clearly, the action of go in 

M^fc naturally extends to M\k- 

Define also the restricted dual module M' x k to M\^ k - the direct sum of the 
dual spaces to the homogeneous subspaces in M\^, with the action of go defined 
by duality. This module is a lowest weight module. Denote its lowest weight vector 
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Now let us define evaluation modules. Let \i e C. Define the module over g 
to be the pullback of the module over gl N defined in Section 1 with respect to 
the evaluation map tt : go — > gl N given by %(A(t) + ac) = A(l) 

Now let us introduce vertex operators. Throughout this paper by a vertex op- 
erator we mean a go- intertwining operator $ : — > Mx,k<S>Vfj,, where Mx^^V^ 
denotes the space of formal sums of the form J2JLo u j ® v j-> u j e Mx,k[—j], Vj £ V^. 

Proposition 3.1. If Mx,k is irreducible then there exists a unique nonzero vertex 
operator up to a scalar factor. 

Proof. In order to construct we first need to construct the vector w = $>vx,k £ 
Mx,k®V^. This vector has to be annihilated by the subalgebra g$ and satisfy the 
condition hiW = \w. To construct such a vector is the same as to construct a 
map w : M' x k — > V M which commutes with the action of an d maps v' x k to a 
zero weight vector. Since M^fe is irreducible, M' x k is freely generated by v' x k over 
C/(0o"), which implies that the map w is uniquely determined once we know the 
image of v' x k under this map. Because the space of zero weight vectors in is 
one-dimensional, w is defined uniquely up to a factor. 

Finally, since $ is an intertwiner, and M\ t k is freely generated by vx,k over 
£7(gQ~), $ is uniquely determined once we know the image of vx,k under it, i.e. once 
we know w. Therefore, $ is unique up to a factor. ■ 

Let A : Mx^k Mx,k be any linear operator. It is obvious that we can uniquely 
represent the operator A as an infinite sum 

(3.4) A=Y,A[m], 

m£Z 

where A[m] = Q)jA[m]j, and A[m]j : Mx, k [j] — > Mx,k\j + m] are linear operators: 
for a homogeneous vector u of degree j in Mx,k the vector A[m] jit is defined as the 
degree j + m component of Au. 

Let Zj, 1 < j < N, and q be formal variables. Let z = (z±, zn)- Then we can 
consider the formal series 

N N 

(3.5) Fx, k (Az,q) = Tr\ Mx , k (Aq~ d J] z)*) = ]T q m Tr \ Mxk[m] (A[0] m l[z^). 

j = l m>0 j = l 

Every coefficient of this formal series is a function of the form Ylf=i z j 3 'Pi^i ■■■■> ~^T~)' 
where p is a Laurent polynomial. 

Now let us define the (modified) expectation value (or 1-point correlation func- 
tion) of an operator A as follows: 

(3.6) < A > x ,k (z) = -=, , , y P = (Pi, Pn), Pj = — J- 

Then < A > x , k e C[(z 1 /z 2 ) ±1 , (z JV _i/z N ) ±1 ] [[q]] for any A. 

Let 9 : V M — > C be a zero weight linear function. This function is unique up 
to a factor. Let the operator $o : ^\,k ~^ ^\,k be defined by the condition: 
= (Id <8> 0)(3>u), u G Mx, k - The main object of our study will be the 1-point 
correlation function of this operator: 
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Finally, let us define the algebra of regular expressions, U(go) veg . 

Let X±(j) be an eigenbasis for {hi} and d in q^, respectively, such that the 
eigenvalue of d on X±(j) is a monotonic function of j. Let us call an expression of 
the form 

(3.8) Y = X_(m) n ™ . . .X_(l) ni /if . . . h p N N X+(l) ri . ..X + (s) r °c M G U(q ) 
a standard monomial. Let us call the number 

m N s 

(3.9) ord(y) =5>i + E^ + E r i + M 

j=i j=i j=i 

the order of the standard monomial Y, and the number 

s 

(3.10) sdeg(Y) = J>,deg(X+) 

i=i 

the subsidiary degree of Y. According to the Poincare-Birkhoff-Witt theorem, the 
standard monomials form a basis in U(qo). 

Definition. 1) A sum Y^m^i^rnYm, where b m G C* and Y m are distinct standard 
monomials, is called a regular expression of order <n if ord(Y m ) < n for all m. 

2) A series Ylm=i ^>mXm, where b m G C* and Y m are distinct standard monomi- 
als, is called a regular expression of order < n if 

(i) ord(Y m ) < n for all m > 1; 

(ii) lim^oo sdeg(Y m ) = oo. 

3) A regular expression is of order n if it is of order < n but not of order < n — 1 . 

It is easy to show that regular expressions form an associative algebra. This 
algebra is a completion of U(go) and is denoted by U(Qo) reg . 

The main property of the algebra U(Qo) reg is that the action of U(qq) in Verma 
modules over $jo can be naturally extended to U(go) reg : if R = Ylm=i ^>mXm is a 
regular expression then for any w G Rw = Ylm=i ^mY m w if n is large enough. 

4. The mapping \ 

Let T>On be the algebra of differential operators in zi, zn with Laurent poly- 
nomial coefficients (over C). Let DO N = VO N [[q]]. Let A = End(V At ) ® DO N . 

We will consider expressions of the form < $Y~ >A,fe) where Y is a linear combi- 
nation of standard monomials. Such expressions are (formal) functions of z, q with 
values in V M . For the sake of brevity, we will write < A > instead of < A >\ : k- 

Theorem 4.1. Let Y be a finite linear combination of standard monomials of order 
< n, and let k be a complex number. Then: 

1) There exists a differential operator Ly(k) G A such that for all A for which 
is irreducible 

(4.1) < $y >= Ly(fc) < $ >, 

satisfying the following conditions: 
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3) The operator Ly(k) is a polynomial in k of degree < n; 

4) If Y is a standard monomial and sdegiY) = s then Ly = q s Ly, where 
L Y e A. 

Proof. The proof of of the theorem is by induction in the order of Y. More precisely, 
we will construct Ly for Y of order n from Lj for standard monomials J of lower 
orders. 

For order 0, Y = 1, and we can take L = 1. Assume that Y is a standard 
monomial of order n > 0. Consider the last factor X e go of this monomial, so that 
Y = Y'X, where Y' is a standard monomial of order n — 1. Since ord(Y') = n — 1, 
we already know that there exists an operator Ly> of order < n — 1 which is a 
polynomial in of degree < n — 1, such that < <&Y' >= Ly/ < $ >. 

There are three possibilities for X. 

1. X = c. In this case we have < &Y >= k < $>Y' >. Then we can set 
Ly = /cLy/. This operator satisfies (4.1). It also satisfies statements 2,3,4 of the 
theorem because so does Ly. 

2. X = X±[j]. In this case let T = [Y'X]. It is clear that T is a finite linear 
combination of standard monomials of order < n — 1, so we may assume that 
the operator Lt is already defined. Then, using the intertwining property of $ 
($X = (X <g) 1 + 1 ® we can write 

(4.2) 

< $y >=< <$>Y'X >=< $XF' > + < $T >=< X<$>Y' > +tt(X) < > + < $T > 
On the other hand, we have 



< X$Y' >= 



(4 3) («^ n ^*) = ft 2? , < $y , . 

where are integers defined by [hi,X] = a{X. Relations (4.2), (4.3) show that 

N 

(4.4) < $Y >= q~ de ^ Y[ < $Y > +(n(X)L Y/ + L T ) < $ >, 

i=l 

which implies that 

(4.5) < $Y >= ^ (Tr(X)Ly/ + L T ) < $ >, 

where denotes E P > ? mp ^ P if m > °> and ~ E P >i 3~ mp £~ p if m < 0. This 

means that the differential operator 

( 4 - 6 ^ Ly = i JmiV a. ^W^ + L T) 



satisfies (4.1). The fact that it satisfies statements 2,3,4 follows from the validity 



13 



3. X = hi. Then we have 



(4.7) < $y >= ^- -§^{F- P ,k < *Y' >), 

which implies that the operator 

(4-8) L Y =-^^-.F_ p , k L YI 

satisfies (4.1). The fact that it also satisfies statements 2,3,4, as before, follows 
from the validity of these statements for Ly . 
The theorem is proved. ■ 

Theorem 4.2. Let Y be a regular expression of order < n and weight (i.e com- 
muting with hi for all i), and let k be a complex number (k ^ Q). Then: 

1) There exists a unique differential operator Dy(k) e DOn such that for all A 
for which is irreducible 

(4.9) < $ Y >= D Y (k) < $ >, 

satisfying the following conditions: 

2) The operator Dy(k) is of order < n as a differential operator; 

3) The operator Dy(k) is a polynomial in k of degree < n. 

4) If Y is a standard monomial and sdegiY) = s then Dy = q s Dy, where 
by e VO N . 

Proof. First assume that Y is a finite linear combination of standard monomials. 
Since Y is of weight 0, the operator Ly preserves weight in V M . Let Dy denote the 
restriction of Ly to the zero weight subspace in V M . Obviously, Dy e DOn, and it 
satisfies properties 2,3,4. In particular, because it satisfies property 4, the definition 
of Dy can be extended to infinite regular expressions: if Y = Ylm=i ^rnXm, where 
Y m are standard monomials, then we can set 



(4.10) d y=J2 b ™ DY ~ 



m=l 

and this series will be convergent as a formal series in q since any fixed power of q 
occurs in only a finite number of its terms. 

It remains to prove that Dy satisfying (4.9) is unique for any Y. Suppose Dy 
and D ( y ] both satisfy (4.9), and let D = D ( y } - D ( y } . Then D < $ >\, k = for 
all A. Assume that D ^ 0. Let D m q m , D m e T>On : be the leading term in the 

^-expansion of D. Since < $ >A,jt|g=o= IljLi ^j J ~ k/N , we have: Andlj^) = 
for generic complex Ai,...,Ajv. This implies that D m = 0, which contradicts the 
assumption that D m is the leading coefficient. Therefore, D = 0, and hence Dy = 

Now we are in a position to define the mapping \. From now on the notation 



14 



First we define the linear mapping \ '■ ^(flo)reg — *• DOn which acts according to 
the formula: 

(4.11) x(Y)=D Y (-N). 

It is easy to see that D c y = kDy , which implies that x kills the ideal (c + 
N)U(go)® eg , and therefore descends to a mapping 

(4.12) X ■■ U(Qo)° reg /(c + N)U(Q )° Teg - DO N . 

This mapping will be the main tool in our proof of the integrability theorem. 

For brevity we denote the algebra U (g ) Ye g/ (c + N)U(q ) reg by U cr i t . This nota- 
tion is suggested by the physical terminology: representations of g with c = — N 
in which U cr it acts are called critical level representations. 

The mapping x has the following important property: 

Lemma 4.3. Let C G U cr {% be a central element: \C,X\ = for any X e U cr n. 
Then 



(4-13) x(CY) = X (Y)x(C), Y e V% 



crit • 



Proof. Let Y, C be representatives of Y,C in U{Qo)^. e& . Then for any X e So 
[C,X] = (c + N)X', X' E U(qq)®. Therefore, it follows from the construction of 
the operators D Y (cf. proofs of Theorems 4.1, 4.2) that 

D d y < $ >=< §CY >= D Y <$C> +(k + N)E < $ >= 

(4.14) D Y D G < $ > +(k + N)E < $ >, Ee DO N , 

i.e. (by the uniqueness of D Y ) 

(4.15) D dY (k) =D Y D d (k) + (k + N)E(k). 
Therefore, 

(4.16) D dY (-N)=D Y D d (-N), 

which implies (4.13), Q.E.D. ■ 
Lemma 4.3 implies 

Theorem 4.4. Let Z be the center ofU cr a. Then the linear mapping x '■ Z — > DOn 
is a homomorphism of algebras, and hence x(Z) is a commutative subalgebra in 
DO N . 
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5. Sugawara operators and the integrability theorem 

It is well known that the algebra U cr [ t has a big center Z. This center is generated 
by the so-called Sugawara operators. Explicit construction of these operators is very 
technical, but for the proof of the integrability theorem we will need a very limited 
amount of information about them. 

Let £1, ...,£jv_i is any basis of the Cartan subalgebra of sIn orthonormal with 
respect to the inner product < A, B >= Tr AB. We assume that £,jt m , m^O, are 
among the basis vectors Xf. Also, if a G sIn, we denote at m as a[m]. 

Sugawara operators for go have orders 1, 2, N. We denote the Sugawara op- 
erator of order j by Qj. The explicit expressions for fli and O2 are very simple: 

N 

Qi = y~]hj, 

i=i 

N 

i=l ijtj 
00 N — l 00 

2 E E 6[-HCiH + 2 ^^[-ml^lm: 

m=l i=l m=l i^j 

About higher order Sugawara operators we need to know the following: 

Theorem 5.1. ([*], Proposition 3.3; see also [*]) There exist elements Tj e Z, 
3 < J ' < iV, 0/ orrfer j, which can be written in the following form: 

N 

(5.2) T j = S j + P j (S lt ...,S j - 1 ) + E j , S m = Y J K 

i=\ 

where P is a polynomial and Ej contains only standard monomials of positive sub- 
sidiary degree. 

Corollary 5.2. There exist degree zero elements Qj<EZ3<j<N, of order j , 
which have the following form: 

(5.3) Q j = S j + E' j , 

where E'- contains only standard monomials of positive subsidiary degree. 
Proof. Define the elements T- 1 ^ of Z as follows: 

(5.4) T^ =T J -P J (T 1 ,...,T J _ 1 ). 
Then we have 

(5.5) T\ l) = S 3 + if } (Si, Sj- 2 ) + Ef\ 

where P- is a polynomial, and all standard monomials in • have positive sub- 
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Next, define Tj 2) G Z by: 

(5.6) Tf } = if } - P,(T 1 (1) , ...,Tji ) 1 ). 
Then we have 

(5.7) T^ 2) = S J +P^(S 1 ,...,S J . 3 )+Ef , 

(21 (21 

where P- is a polynomial, and all standard monomials in E y - ' have positive sub- 
sidiary degree. 

Continuing this procedure, we will get to the j — 1-th step, and define T^~ X \ 

which will have the decomposition Sj + E^~ X \ Therefore, we can set fij = T^~ X \ 
Q.E.D. ■ 



Let us now assume that Zj, q are complex numbers such that Zj ^ 0, < \q\ < 1, 
and introduce new variables Xj G C and r G C + such that e 27riXj = Zj, e 27riT / N = q. 
For the proof of the integrability theorem we need one more result: 

Theorem 5.3. 

(5.8) X(^ 2 ) = + const ' 
where H is defined by the formula 

(5.9) ^ = ^ 2^ + 1) Pi^-Xj-^rlr). 

j=l j l<i<j<N 

Proof. It is shown in [EK, Proposition 4.2] that ip\ t k satisfies the following parabolic 
differential equation: 

(5.10) -47v 2 (k + N)q^- = (H + c)^ k , 

where c is a constant. But we know that (k + N)q^^- =< <& (k + N)d >=< 

$0^2 >= Dq 2 < $ >= Dn 2 ip\,k, which implies that Dq 2 = (An 2 )' 1 !! + const, 
Q.E.D. 

Now we are ready to prove the complete integrability theorem for the elliptic 
case. 

Theorem 5.4. f [OP],) The Hamiltonian H given by (1) with U = Kp is completely 
integrable. 

Proof. . Let fx be such that K = — 2//(rmt + 1). Consider the differential operators 
Df ll ,...,Dci N . They are pairwise commutative, and, thanks to Corollary 5.2, 

■A & 

(5.11) Dq, = (2-7T c) J > — h lower order terms 

i=i 3 

Let us make a change of variable Xj = Xj + j^r, and let Lj be the image of (2% C 
) J Dq. under this change of variable. Then [Lj, L m ] = for all j, to, L 2 = H + const 

because of Theorem 3.3, and Lj = (m^) ^ l° wer order terms, which implies 
that $Tja\ rtp a.lfffihra.ir.a.llv inHfvnfmdfmt. Thp thmrfim is nrnvfiH. 
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